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Data from e + e — annihilation into hadrons collected by the JADE experiment at centre-of-mass en- 
ergies between 14 and 44 GeV were used to study the 4-jet rate using the Durham algorithm as well 
as the first five moments of event shape observables. The data were compared with NLO QCD pre- 
dictions, augmented by resummed NLLA calculations for the 4-jet rate, in order to extract values of 
the strong coupling constant ag. The preliminary results are ag(Af z o) = 0.1169 ± 0.0026 (4-jet rate) 
and ag(Af z o) = 0.1286 ± 0.0072 (moments) consistent with the world average value. For some of the 
higher moments systematic deficiencies of the QCD predictions are observed. 



1 Introduction 

The production of hadrons in e + e _ annihila- 
tion allows precise tests of the gauge theory 
of strong interactions, Quantum Chromody- 
namics (QCD). In this paper recent and pre- 
liminary analyses of JADE data using the 4- 
jet rate based on the Durham algorithm [1] 
and using the first five moments of event 
shape observables are presented [2,3]. 

The data used in our analyses were col- 
lected at centre-of-mass (cms) energies -y/s = 
14.0, 22.0, 34.6, 35.0, 38.3 and 43.8 GeV be- 
tween 1981 and 1986 with the JADE detec- 
tor [4]. The data samples consist of 0(1000) 
events at yi = 14.0, 22.0, 38.3 and 43.8 GeV 
while at Vs = 34.6 (35.0) GeV about 14000 
(21000) events are used. 

The software employed to perform the 
analyses includes the original JADE detector 
simulation, event reconstruction and event 
display programs, see [2,3] for details. It is 
possible to use recent Monte Carlo event gen- 
erators such as PYTHIA, HERWIG or ARI- 
ADNE to generate simulated events, to pass 
these through the JADE detector simulation 
and to reconstruct them in essentially the 
same way as the data. The event generators 
were used with parameter settings obtained 
by OPAL after adjusting to LEP 1 data. The 
original data are only available after a further 
step of data reduction has been performed, 



resulting in information about 4- vectors of re- 
constructed particles and some quantities for 
event selection. 

The selection of well reconstructed tracks 
from the tracking detectors and clusters from 
the electromagnetic calorimeter as well as 
the selection of hadronic events follows previ- 
ous analyses. Most importantly requirements 
on visible energy and momentum, balance 
of momentum along the beam direction and 
track multiplicity suppress events from two- 
photon interactions, r production and other 
backgrounds [4-7]. 

The contribution of e + e~ — * bb events 
is subtracted from the data using simulated 
events. Then the data are corrected for the 
effects of detector resolution and acceptance 
of selection cuts by correction factors deter- 
mined from simulated events before and after 
the JADE detector simulation. 

Experimental systematic uncertainties of 
the analyses include variation of the event 
selection cuts, comparison of data sets re- 
sulting from different versions of the JADE 
event reconstruction program and comparing 
experimental corrections derived from differ- 
ent event generators. 

Additional systematic uncertainties are 
studied for the extraction of values of the 
strong coupling constant as- Hadronisation 
uncertainties are evaluated by using differ- 
ent Monte Carlo generators to compute the 
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hadronisation corrections. Theoretical sys- 
tematic uncertainties are found by changing 
the renormalisation scale \i of the QCD pre- 
dictions from fj, = y/s to fj, = y/s/2 and 

fj, = 2y/s. 

2 4-Jet Rate 

Jets are reconstructed in the hadronic 
events using the Durham jet clustering algo- 
rithm [1]. Figure 1 shows the data for the 4- 
jet rate corrected for experimental effects as 
a function of y cut measured at y/s = 35 GeV. 
The data are compared to predictions from 
the event generators PYTHIA, HERWIG and 
ARIADNE. We find good agreement between 
the data and the predictions within the sta- 
tistical and experimental errors and conclude 
that the models can be used to correct for 
hadronisation effects. 
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Figure 1. The figure shows the corrected data for the 
4-jet rate at \fs = 35 GeV compared with predictions 
from Monte Carlo models. The error bars represent 
statistical and experimental errors added in quadra- 
ture [2]. 



The QCD predictions arc next-to- 
leading-order (NLO), i.e. 0(a|) in leading- 
order (LO) with 0(a.g) radiative corrections, 
combined with resummed next-to-leading- 
logarithm (NLLA) calculations [8]. We find 
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Figure 2. The results for as from the 4-jet rate are 
shown as a function of y/s. The error bars show the 
statistical and total uncertainties. The full and dash- 
dotted lines indicate the current world average value 
of a s (M z o ) [9] . The results at yfs = 34.6 and 35 GeV 
have been combined for clarity. Results from OPAL 
and ALEPH are shown as well [2]. 



good agreement between data and theory for 
large values of y cut = 0(1O~ 2 ) where mostly 
3- and 4-jet events are found and the theory 
considering at most five partons is expected 
to be reliable. The resulting values of as are 
shown in figure 2. 

The values for as (Af z o ) from the data at 
s/s = 22.0, 34.6, 35.0, 38.3 and 43.8 GeV are 
combined taking into account correlations be- 
tween experimental, hadronisation and the- 
oretical systematic uncertainties. The fits 
at y/s = 14.0 GeV have large experimen- 
tal and hadronisation uncertainties and are 
therefore excluded from the average. The re- 
sult is a s (M z o) = 0.1169 ± 0.0004(stat.) ± 
0.0012(cxp.) ± 0.0021(had.) ± 0.0007(thco.), 
a s (M z o) = 0.1169 ± 0.0026 (total error), 
consistent with the current world average 
a s (M z o) = 0.1182 ± 0.0027 [9]. 
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3 Moments of Event Shape 
Observables 

The first five moments of the distributions of 
the event shape observables 1 — T, C, Bt, 
Bw, J/23 and Mh are calculated according to 
(y n ) = Jo'™'" 1 V n l/^dcr/dj/dy', where y de- 
notes one of the observables, y max is the kine- 
matically allowed upper limit of the observ- 
able and n = 1, . . . , 5. 

The calculation of perturbative QCD 
predictions in NLO (£>(a s ) in LO with 0{a%) 
radiative corrections) involve a full integra- 
tion over phase space. This analysis is thus 
complementary to tests of the theory using 
the differential distributions which are com- 
monly only compared with data in restricted 
regions, where the theory is able to describe 
the data well, see e.g. [5]. 

Figure 3 presents the data for the first 
five moments of 1 — T, Mh and C corrected 
for experimental effects compared with pre- 
dictions by the same event generators as in 
section 2. There is generally good agreement 
between data and model predictions; HER- 
WIG is seen to describe the data somewhat 
less well than PYHTIA or ARIADNE. We 
will use the models to derive hadronisation 
corrections in order to compare the data with 
predictions from perturbative QCD. 

We fitted the QCD predictions corrected 
for hadronisation to the data for a given ob- 
servable and moment n = 1 , . . . , 5 individu- 
ally with as(M z o) as the only free parame- 
ter. The results for as(M z o) are summarised 
in figure 4. The fit to (Mh) did not con- 
verge and therefore no result is shown. We 
observe that the values of as(M z o) increase 
with n for the observables ((1 - T) n ), (C n ) 
and (BJ}.), while for the other observables 
(#w), (G/23)") and (Mg), n = 2,..., 5, the 
results are fairly stable. 

We evaluated the ratio K of NLO and 
LO coefficients for the six observables used in 
our fits and found a clear correlation between 
the steeply increasing values of as (M z o ) and 
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Figure 3. The corrected data for the first five mo- 
ments of the observables 1 — T, Mjj and C are pre- 
sented with error bars showing statistical and experi- 
mental error added in quadrature. The lines indicate 
the predictions of Monte Carlo models. The lower 
panels show the differences between data and models 
at y/s = 14 and 35 GeV, divided by the errors [3] . 



increasing values of K with n for ((1 — T) n ), 
(C n ) and (B%). The other observables (B^), 
((2/23)™) and (M™), n = 2, . . . , 5, have fairly 
constant values of K and correspondingly 
stable results for a s (M z o). We also noted 
that (Mh) has a large and negative value of 
K which is the cause that the fit did not con- 
verge. 

In order to find a combined value of 
as (M z o ) we considered only those results for 
which the NLO term is less than half the LO 
term (i.e. \Koi$/2-k\ < 0.5), namely (1-T), 
(C), (B T ), TO and ((t/23)"), n = 1.....5 
and (Mg), n = 2, . . . , 5; i.e. results from 17 
observables in total. The purpose of this 
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Figure 4. Measurements of ag(M z o) using fits to mo- 
ments of six event shape observables are shown. The 
inner error bars represent statistical errors, the mid- 
dle error bars include experimental errors and the 
outer error bars show the total errors. The dotted 
line indicates the weighted average described in the 
text; only the measurements indicated by solid sym- 
bols were used for this purpose [3]. 



requirement was to select observables with 
an apparently converging perturbative pre- 
diction. Correlations between statistical, ex- 
perimental, hadronisation and theoretical un- 
certainties were considered when forming the 
average. The result is a s (M z o) = 0.1286 ± 
0.0007(stat.) ± 0.0011(exp.) ± 0.0022(had.) ± 
0.0068(theo.), a s (M z o) = 0.1286 ± 0.0072 
(total error), above but still consistent with 
the world average value. It has been observed 
previously in comparisons of distributions of 
event shape observables with NLO QCD pre- 
dictions with rcnormalisation scale /i = ^/s 
that fitted values of as(M z o) tend to be large, 
see e.g. [10]. 



4 Summary 

We have presented preliminary results of 
measurements of the 4-jet rate based on the 
Durham algorithm and the first five moments 
of event shape observables using JADE data 



at s/s = 14.0 to 43.8 GeV. The predictions of 
the Monte Carlo models PYTHIA, HERWIG 
and ARIADNE tuned by OPAL to LEP 1 
data were found to be in reasonable agree- 
ment with the data. The data have also been 
used to extract measurements of the strong 
coupling constant as(M z o) with the results 
a s (M z o) = 0.1169 ± 0.0026 (4-jet rate) and 
a s (M z o) = 0.1286 ± 0.0072 (moments) in 
agreement with the world average value. The 
higher moments of 1 — T, C and £?t are ob- 
served to yield systematically larger values of 
a s (M z o). 
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